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Abstract 

Flaton models of Peccei-Quinn symmetry have good particle physics motivation, 
and are likely to cause thermal inflation leading to a well-defined cosmology. They can 
solve the /i problem, and generate viable neutrino masses. Canonical flaton models 
predict an axion decay constant FpQ ~ 10^'^ GeV and generic flaton models give FpQ > 
10^ GeV as required by observation. The axion is a good candidate for cold dark matter 
in all cases, because its density is diluted by flaton decay if -Fpq > 10^^ GeV. In addition 
to the dark matter axions, a population of relativistic axions is produced by flaton 
decay, which at nucleosynthesis is equivalent to some number 6Ni, of extra neutrino 
species. Focussing on the canonical model, containing three flaton particles and two 
flatinos, we evaluate all of the flaton-flatino-axion interactions and the corresponding 
axionic decay rates. They are compared with the dominant hadronic decay rates, 
for both DFSZ and KSVZ models. These formulas provide the basis for a precise 
calculation of the equivalent dNi^ in terms of the parameters (masses and couplings). 
The KSVZ case is probably already ruled out by the existing bound 5N^ < 1. The 
DFSZ case is allowed in a significant region of parameter space, and will provide a 
possible explanation for any future detection of nonzero 5N^. 



1 Introduction 



With the discovery of the instantons it was reahzed that the ^qcd parameter of the Standard 
Model can have important physical consequences. In particular, the induced CP violation 
affects the electric dipole moment of the neutron, leading to the upper limit limit ^qcd < 
10~^°. An attractive explanation for such a small parameter is provided by the Peccei-Quinn 
mechanism [|T], |^. There is supposed to be a spontaneously broken global U{1) symmetry 
(PQ symmetry), which is also explicitly broken by the color anomaly. Its pseudo-Goldstone 
boson is the axion. The PQ symmetry is spontaneously broken by some set of scalar fields 
(pi (elementary or composite) with charges Qi, so that their PQ transformation is 

0. ^ e^Q."0^ . (1) 

Denoting the vacuum expectation value by fj/v^, we define the PQ symmetry breaking 
scale FpQ by 

^^Q = EQk^ (2) 

i 

(In defining -FpQ, we use the canonical normalization of the PQ charges, that the smallest 
is set equal to 1.) Collider and astrophysics constraints require roughly |^ 

FpQ/iV> lO^GeV. (3) 

The bound is actually one on the axion mass, whose relation to Fpq is given by 

ma = 6iV X 10-^ eV . (4) 

Here, is the number of distinct vacua, or equivalently the number of domain walls meeting 
at each PQ string. 

PQ charge will also be carried by fields which do not spontaneously break PQ symmetry. 
We shall consider KSVZ (hadronic) models in which these are only some extra heavy quark 
superfields, and DFSZ models in which they are only Standard Model (SM) superfields. 

We are concerned[| with models in which the fields breaking PQ symmetry are flatons. 
Flatons are fields whose tree-level potential is flat in the limit of unbroken renormalizable 
supersymmetry, which acquire nonzero vevs after soft supersymmetry breaking. We make 
the usual assumption, that the supersymmetry breaking is gravity-mediated. 

In the rest of this section we recall the essential features of fiaton and non-flaton models 
of PQ symmetry. In the next section we discuss in some detail the cosmology of fiaton 
models. In Section |^, we give the general structure of the fiaton and flatino masses. In 
Section ^ we analyze the general self interactions between flatons and fiatinos. In Section 
^ we see the effect of the interaction of the flatons with the matter fields. In Section ^ we 
consider the implication of the nucleosynthesis bound on the energy density of relativistic 
axions, in particular, of a possible future bound SN^ < 0.1. We conclude in Section |^. 

^An earlier version of this work appeared as |iep-ph/9903286. 
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1.1 Models of PQ symmetry-breaking 



Let us consider the potential of the fields (pi which break PQ symmetry. In a supersymmetric 
model there have to be at least two, but for first orientation we pretend that there is only 
one. In the limit of exact PQ symmetry, its potential will be of the form 

V = V,-m'W + -A|0r + E A.^^ • (5) 

n=l ^'^Pl 

Mpi = {MGnY^''^ = 2.4 X lO^s GeV is the reduced Planck mass. 



1.1.1 Non-fiaton models 

If the renormalizable coupling A is of order 1, the non-renormalizable terms are negligible 
and the mass of the radial oscillation is m ~ -FpQ- (Recall that FpQ = ^/2Q {\(f)\) .) 

Going on to the case where a number of fields (pi {i = 1, ■ ■ ■ ,p) break PQ symmetry, there 
will be p particles corresponding to the radial oscillations, and (in addition to the axion) 
p — 1 particles corresponding to the angular oscillations. There will also be p superpartners 
with spin 1/2. In a non-supersymmetric theory all of these particles would have mass of 
order -Fpq, but supersymmetry protects the mass of one scalar and one spin-1/2 particle. 

Indeed, in the limit of unbroken supersymmetry, the holomorphy of the superpotential 
ensures that PQ symmetry is accompanied by a symmetry acting on the radial parts of the 
PQ charged fields 

0. e«'>, . (6) 

The corresponding pseudo-Goldstone boson is called the saxion (or saxino), and the spin- 
1/2 partner of the axion/saxion is called the axino. With gravity-mediated supersymmetry 
breaking, the saxion and axino will typically both have soft masses of order 100 GeV, al- 
though specific models exist with an axino mass of order keV. In gauge-mediated models, 
the saxion and axino will typically both have the sub- keV mass. 

In non-fiaton models one can hope to understand a value FpQ ~ 10^'^ GeV since that is 
the supersymmetry breaking scale 0, but it may be hard to understand a bigger value. 



1.1.2 Flaton models 

We are concerned with models p], ^ |10|, |ll|, in which the fields breaking PQ sym- 
metry are fiaton fields 
of unbroken 



This means that their tree-level potential is fiat in the limit 
renormalizable supersymmetry. We assume that supersymmetry breaking is 
gravity-mediated, which is usual in the context of fiaton fields. 

Considering first the case of one fiaton with potential Eq. (^, the quartic term is (prac- 
tically) zero, while the mass m of the flaton is of order 100 GeV to 1 TeV. When making 
numerical estimates, we shall take for definiteness 



m 



10 



2.5±0.5 



GeV. 



(7) 
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The leading non-renormalizable term generates a large vacuum expectation value (vev), 

FpQ = V2{\<P\) ~ [mM^J^^ . (8) 

Here, A„ is the coefficient of the leading term in Eq. which is expected to be roughly of 
order 1. In the case n = 1 this gives a vev of order 10^° GeV, but it can be bigger if n is 
bigger. Imposing the condition that V (practically) vanishes in the vacuum gives the height 
of the potential, 

1^106 GeV ) [W^GeVj ' ^ ^ 

We are concerned with the case that the (pi are fiaton fields. The scalar particles are now 
the axion, plus 2p — 1 flatons with mass of order 100 GeV, and p flatinos with masses of the 
same order. The saxion, defined through Eq. (^), is a linear combination of flatons, while 
the axino (defined as the partner of the axion-plus-saxion) is a linear combination of flatino 
mass eigenstates. Neither of them has any special significance. In particular, the possibility 
of a keV-mass axino does not exist in fiaton models. In the models that we shall consider, 
p = 2 so that there are three flatons and two flatinos. 



1.2 Estimates of Fpq 

To estimate the magnitude of Fpq in fiaton models, we begin with schematic case of one 
fiaton field, with potential Eq. (||). 

Consider first the expected magnitude of the coefficients A„ in Eq. (|^). In the case 
that Mpi represents the ultra-violet cutoff for the effective field theory, one usually assumes 
A„ ~ 1, but A„ ~ l/(2?T, + 4)! may be more realistic [ITl]. On the other hand, one might quite 



reasonably suppose that the cutoff, call it Auv, is around the gauge coupling unification 
scale 10~^Mpi (either because fields of a Grand Unified Theory have been integrated out, or 
because this is the true quantum gravity scale). In that case, the estimate A„ ~ l/(2n + 4)! 
to 1 would be reasonable if Mpi were changed to Auv in Eq. (H). Retaining Mpi, one should 
multiply the estimate of A„, by a factor (Mpi/Auv)^" ~ 10'^". In view of these considerations, 

we adopt as a reference the estimate l/(2n + 4)! < A„ < 10"^", corresponding to Af = io°-^='=°-^ 

and A| = loO-3±i-o. 

Using these estimates of A„, we can make estimates of Fpq bearing in mind the uncer- 
tainty Eq. (1^) in m. In these and other estimates, we add in quadrature different uncertainties 
in the exponents. This procedure has no particular basis, but at least it is better than ignor- 
ing the uncertainties completely, or adding different estimates linearly. In the case at hand, 
the uncertainty is dominated by the large uncertainty that we assigned to the A„. We take 
the PQ charge of to be 1, so that Fpq = ^/2{\(f)\). 

Unless it is forbidden by a symmetry, the leading term n = 1 will be the one appearing 
in Eq. (§), leading to 



FpQ = 10^0A±0.9 Qgy (n = 1) . (10) 
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If the leading term is n = 2 this becomes FpQ = iQ^'^-^->^^-^ GeV. At higher n, FpQ slowly 
increases, so for n > 1 we have 

FpQ > 10^^-^ GeV (n>l). (11) 

In the class of supersymmetric models that we shall consider, PQ symmetry is broken 
by two fiaton fields 0p and 0q, with charges respectively 1 and 2n + 1. The fields interact, 
giving the rather complicated potential Eq. (RO) for n = 1 and one of similar form |]TB[ for 



bigger n. In all cases, the vevs of both fiaton fields are given roughly by Eq. @. There are 
additional uncertainties because there are more soft parameters and factors involving n, but 
in view of the large uncertainty we already assigned to the coupling A„ Eqs. (|10D and (|Tl]) 
should still provide reasonable estimates. 



2 Cosmology and dark matter 



2.1 Cosmology of the PQ fields 

Generic models of PQ symmetry have many possible cosmological consequences, which have 
been discovered gradually over the years. 

In all cases the axion lifetime is longer than the age of the Universe p[ , so that its present 
density must be f2a ^ 0.3 or so, with the equality prevailing if axions are the dark matter. 
The density depends on whether the axions come from strings or from the vacuum fluctuation 
of the axion field during inflation. (In the latter case the axion density [16| depends on our 
location in the universe.) Strings may be produced by a variety of mechanisms during |1^ 
and after inflation. The axion density depends also on the amount of any entropy production 
after the axion mass switches on at a temperature around 1 GeV. As a result there is no 
model-independent prediction for f2a as a function of Fpq. 

Within the usual framework of non-flaton models, the superpartners of the axion can 
also have a range of cosmological consequences. An axino with a keV mass is a dark matter 
candidate, which may be produced by a variety of mechanisms and give rise to a variety of 
cosmological consequences [jl8|, Alternatively, an axino with a 10 GeV mass may be the 
cold dark matter [|l9l, as it can be the lightest supersymmetric particle (LSP). Finally, the 
saxion is a late-decaying particle which may be produced by thermal or other mechanisms. 
If it is sufficiently abundant to dominate the density of the Universe, it must decay well 
before nucleosynthesis, and before it does so it will dilute the abundance pre-existing relics, 
including baryons and dark matter candidates pO|, If it is less abundant it may decay 



much later, and affect the formation of large-scale structure p2 



In contrast with this generic situation, the cosmology of fiaton models is rather well- 
defined, on the reasonable assumption PBI, [TH] that the PQ fiaton fields generate an era 



of thermal infiation |2|, 25, 23, 26, |2^, 2^, which is not followed by any other such era. 
Thermal infiation occurs long after the ordinary infiation which is supposed to be origin of 
structure, and may wipe out all previously existing relics. When it ends, PQ strings are 
produced, and fiaton decay (the analogue of saxion decay) produces a calculable amount of 
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entropy, with the reheating at a calculable temperature in the range MeV to TeV. As a 
result, the axion density is in principle calculable, and appears to be compatible with the 
observed dark matter density for any FpQ allowed by other considerations (Section p.3| .) In 
other words, the axion is a good dark matter candidate in flaton models. Finally, the flatinos 
(the generalizations of the axino) cannot have the keV mass, and for simplicity we assume 
that none of them is the LSP. (The opposite case will be explored in a future paper [^.) 

A unique feature of flaton models is that fiaton decay creates a highly relativistic popula- 
tion of axions p5| , p3[ . This population has nothing to do with the dark matter. Its density 
at nucleosynthesis is equivalent to roughly SN^, ~ 1 extra neutrino species Flaton mod- 
els of PQ symmetry will therefore be a candidate for explaining a nonzero 6N1, that may be 
established in the future, and the models will be strongly constrained if the present bound 



[301 6N,y ^ 1 is significantly tightened. 



2.2 Thermal inflation and reheating 

In the early Universe, with Hubble parameter H > 100 GeV, fields with the true soft mass 
100 GeV are expected ||3l] to have an effective mass-squared rn?{t) ~ =bif^. (During 



\m 



inflation this result might be avoided |Q, but it should still hold afterwards |34 



This applies in particular to the fiaton fields. Pretending for the moment that there is only 
one fiaton field, we focus on the case that m?{t) is positive in the early Universe, because 



thermal inflation [p4| , |25| , p^ , p6| , P7| , ^ then occurs, leading to rather definite predictions 
for the cosmology. 

Let us summarize the history. After H falls below |m|, there will be enough thermalization 
to hold (j) at the origin until thermal inflation begins. (See Section IIIB of Thermal 
inflation begins when the potential Vq dominates the energy density. This is at the epoch 
T ~ Vq'^ ~ 10^ GeV (assuming for simplicity that full reheating has occurred by that time). 
Thermal inflation ends after ~ In (^Vj/'^'^/lml^ < 10 e- folds, when the temperature is of order 
\m\ ~ 100 GeV. 

When thermal inflation ends, the flaton field moves away from the origin. Cosmic 
strings form, and between them the roughly homogeneous flaton field starts to oscillate 
around its vev. Corresponding to the oscillation is a population of flatons. We discuss 
in Section p.4| the possibility that parametric resonance rapidly drains away the energy of 
this oscillation, finding that this phenomenon will probably not occur and will in any case 
have little effect on the following considerations. Discounting parametric resonance, energy 
loss comes at first only from the Hubble drag, which is negligible during one oscillation. 
The oscillation corresponds to flatons, with conserved number and non-relativistic random 
motion (matter as opposed to radiation). When the flatons decay, the Universe thermalizes, 
at the reheat temperature p3l 



Trh ^ 10°-=^(7RiVMpir~10-°-YlO-2MpicAr,hanm3Fp-c^ (12) 

^GeV\ / m \ 
^PQ I V 102-5 GeV j 
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In the first expression, F is a typical fiaton decay rate, while (^rh ~ 100 is the effective 
number of particle species at Trh- In the second expression A'chan is the number of decay 
channels, c is a factor of order 1 and m is a typical soft parameter. The factor 10~^ in the 
estimate of T is what one expects in the case of unsupressed couplings for a single decay 
channel ^3\, and it is confirmed by, for instance, the estimates in Eq. (|5^) , etc. There are 
in reality several channels, and for definiteness, we take cNchan = 10^ °^^'°, leading to the 
estimate in the third line. 

Adding in quadrature the uncertainty of Eq. (^, we find 



RH 



10--fl^^l GeV. 



(14) 



In order not to upset nucleosynthesis, it is required to have Trh ^10 MeV, which corresponds 
to 

i^PQ< lO^'GeV. (15) 



Using Eqs. (|TOD, ([TT]), (g) and (0), we estimate 



= lO^-^^^-^' GeV 
Trh < 10°-^ GeV 



{n = l] 
(n>l). 



(16) 
(17) 



In this discussion of thermal inflation and its aftermath, we have retained the pretense 
that there is just one fiaton field. There are in the models we shall consider two fiaton fields 
(j)p and (pQ, with the potential Eq. (|30D or its n > 1 analogue. We assume that mp(t) is 
positive in the early Universe, so that thermal inflation occurs. When thermal inflation ends, 
0P moves away from the origin, and as a result 0q also moves away from the origin. At first 
the orbit in field space will be far from the vev, but after a few Hubble times the Hubble 
drag will allow the vev to attract the orbit towards it, so that there are almost sinusoidal 
oscillations of the eigenmodes around their vacuum values. These oscillations are equivalent 
to the presence of the three species of fiaton, and each of them decays at the epoch specified 
by Eq. i\L2\) , with m the appropriate mass. The reheating process is complete after the last 
decay has taken place. 



2.3 Dark Matter and baryons 

Axionic dark matter 

The axion number density is conserved after some epoch Tcons ~ 1 GeV. In the case n = 1, 
Trh is bigger than T^ons and entropy is conserved too. The axion density is then expected 
to be of the form 

V 1012 GeV. 

The constant C is in principle calculable from the dynamics of the strings, walls and axions, 
derived ultimately from the field equation of the fiaton fields breaking PQ symmetry. Ac- 



rr X 1.2 



cording to one group [36, 37, 38], C ~ 1 to 10, while according to another [R9, HO, 41, 42 
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C ~ 0.2. In view of this uncertainty, which we emphasize is one of computation rather than 
principle, we conclude that the axion is a good dark matter candidate in the case n = 1 which 
corresponds to Eq. ([T0|) . (By a 'good' candidate, we mean one whose density is predicted to 
be within at least a few orders of magnitude of the observed dark matter density.) 

In the case n > 1, Trh is smaller than Tcons- Entropy is produced until the epoch Trh, 

0.44 



givmg 




(19) 



with roughly C ~ IOC. At least if C is not too big, the axion is a good dark matter candidate 
in these models too |[15||. 



Baryogenesis and the LSP 

If thermal inflation wipes out pre-existing relics, baryogenesis has to occur after thermal 
inflation. The crucial factor here is the final reheat temperature Trh- Baryogenesis mech- 
anisms occurring at the electroweak phase transition can operate if Trh ^ 100 GeV. This 
is possible in the canonical case n = 1, but not in the case n > 1. If Trh is smaller one 
must turn to other mechanisms, which are quite speculative. Proposals include a compli- 
cated Affleck-Dine mechanism along the lines of QCD baryogenesis |44] or parametric 



resonance baryogenesis [45 



If the lightest supersymmetric particle (LSP) is stable, it has to thermalize in order to 
avoid overproduction from flaton decay. This requires Trh ^ "^lsp/20 [|15|, say Trh more 
than a few GeV. This is exactly what one expects in the case n = 1. As is well known, a 
stable LSP is a good dark matter candidate. 

If the LSP is unstable (due to i?-parity violation), baryogenesis can occur simply by 
allowing a baryon-number violating flaton decay channel (the Dimopoulos-Hall mechanism 
[^ ). This mechanism requires DSFZ as opposed to KSVZ coupling to matter, and as we 
shall see the former case is favored in flaton models. For the mechanism to work, flnal reheat 
must occur at a temperature less than a few GeV This is likely for n > I, but looks 
rather unlikely in the canonical case n = I. 

Let us summarize. In the canonical case n = I, the LSP can thermalize, and therefore 
can be stable so that it is a good dark matter candidate just like the axion. In this case, 
baryogenesis mechanisms involving the electroweak phase transition can operate. In the case 
n > I, the LSP cannot thermalize and therefore cannot be stable. Baryogenesis from flaton 
decay (the Dimopoulos-Hall mechanism) is a natural possibility in this case. 



Supermassive dark matter 

We have seen that the axion is a good dark matter candidate, and that in the canonical 
model the LSP is also a good candidate. These conclusions hold both in the DFSZ and 
KSVZ cases. In the KSVZ case, there is a third good dark matter, namely the heavy quarks 
E,E'^, which are strongly interacting massive particles (SIMPs). Until thermal inflation 
ends, the SIMPs are light and will be in thermal equilibrium. After thermal inflation ends. 
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SIMPs acquire mass through the couphng (j)EE^ to the flaton. The density of such particles 



has been shown p8| to be naturally in the right ballpark. 

The scenario of should be contrasted with an earlier proposal |4^, that the super- 
massive particle is very heavy also in the early Universe, and never in thermal equilibrium. 



Such a particle may be produced by the vacuum fluctuation during ordinary inflation 



or by other mechanisms. In the former case, each comoving wavenumber leaving the horizon 
during inflation will acquire very roughly one particle per quantum state [^8|, |5^ , leading 
very roughly to [^, 



1014 GeV/ VlQi^GeV 



^INFRH \ 

109 GeVy ^' 



(20) 



In this expression, m is the mass of the superheavy particle, H^, is the Hubble parameter 
during slow-roll inflation, and Tinfrh is the temperature at reheat after slow-roll inflation, 
and 7 is the dilution caused by entropy production after that epoch. In our case, thermal 
inflation will give roughly 7 ~ 6^^°. One can adjust the other parameters to make fio ~ 1, 
but in contrast with the LSP and the axion the required value f2o ~ 1 is not favored over 
any other. In other words, this kind of supermassive dark candidate is not (in our present 
state of knowledge) a good dark matter candidate, merely a possible one.|^ 



2.4 Parametric resonance? 

To check whether parametric resonance |5ll| occurs, we make the very crude approximation 



that the sinusoidal oscillation corresponding to the three flatons is present from the very 
beginning. We also assume that the masses of the three flatons are roughly the same, or else 
that one of the amplitudes is much bigger than the others. Then the field equation of the 
Fourier component of each produced field 0„ is a Mathieu equation, leading to a situation 
that has been analyzed in the literature [^. (More realistic cases, including the one where 
the oscillation starts at a maximum of the potential seem to give similar results.) The 



oscillation of the real field 0/ corresponding to the /th flaton leads to an oscillating mass- 
squared m^(0/) for each of the produced scalar fields. Parametric resonance occurs, leading 
to significant production of 0„, if |^ 



g = m^„(0,)/m^ > 10^ (21) 

In this formula, m ~ 100 GeV is a typical oscillation frequency, and mo„(0/) is the amplitude 
of m'^{(j)j). The initial amplitude of oscillation is 0/o ~ -^pq- The scalar particles that can be 
produced include the flatons and the axion, and in the DFSZ case also the Standard Model 
Higgs and sfermions. From Sections § and ^ each of these has mon ~ 100 GeV, making 
g ~ 1. The conclusion is that parametric resonance probably does not occur, but a more 

■^The same is true in the case 7 = 1. In particular, the choice m ~ advocated in is not in fact 
particular favored, bearing in mind that m is the true mass as opposed to the effective mass during inflation. 
Inflation with a potential V = 5"t-inf'?^inf indeed ends when H ~ rrii^p, but the inflaton field (pwp is not 
supposed to be stable and hence is not a dark matter candidate. 
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detailed calculation is needed to say anything definite especially in view of the extremely 
anharmonic nature of the initial oscillation. 

Even if it occurs, the overall effect of parametric resonance will not be dramatic, unless 
it leads to baryogenesis P5|. Its initial effect is to quickly damp the flaton field oscillation. 



producing flatons, axions and in the DFSZ case Higgs and sfermions. The produced particles 
have very roughly the same energy density, and are marginally relativistic except for the 
axions which are highly relativistic. The flatons are stable on the Hubble timescale while 
the Higgs and sfermions decay into highly relativistic ordinary matter (plus the marginally 
relativistic LSP if it is stable). After a small number of Hubble times the energy in the highly 
relativistic particles becomes negligible. The dominant energy is in non-relativistic flatons, 
coming partly from the parametric resonance and partly from the residual homogeneous 
oscillation of the fiaton fields. Except for the baryogenesis possibility, there is no change. 



2.5 Relativistic axions 

Now we come to the relativistic axions, which will be our concern for the rest of the paper. 
This axion population comes from the decay of the flatons H when they finally reheat 
the Universe. Its density during nucleosynthesis is conveniently specified by the equivalent 
number of extra neutrino species 

, (22) 



NS 



where pa is the energy density of relativistic axions, and is the energy density of a single 
species of relativistic neutrino. 

At present, constraints coming from nucleosynthesis are bedeviled by the fact that there 
are two separate allowed regions of parameter space, corresponding to 'low' and 'high' deu- 



terium densities. In the 'high' region, the bound |^0[ at something like 2-a level is AA^,^ < 1.8, 
while in the perhaps favored 'low' region, the bound at a similar level is AA^j^ < 0.3. As we 
shall see, the flaton models predict ^A";^ roughly of order 1, so that at least the second bound 
is quite constraining. 

In the canonical model that we shall discuss, there are three flaton species, and in more 
general models there are more flatons. In general, each flaton species can decay into rel- 
ativistic hadronic matter X, into Xa where a denotes a relativistic axion, or into aa. (We 
neglect for simplicity the small branching ratio into channels containing more axions.) Let 
us pretend first that there is one flaton 0. The hadrons X thermalize immediately, but the 



axions do not thermalize [T^, so their density at reheating is 

Pa = -BaPr , (23) 

where 

^ r(0 ^ aa) + |r(I ^ Xa) 

r(0^x) ' ^ ^ 



9 



and the density of thermalized radiation is 



After reheating, pa oc (ajm/a)'^ = (sns/srh)^, where a is the scale factor, and s is the entropy 
density of the particles in thermal equilibrium. The latter is given by s = {27r'^ / 45) gT^, 
where g is the effective number of relativistic species in thermal equilibrium and T is their 
temperature. At the beginning of the nucleosynthesis era, g = 10.75 and 

TT^ /10.75\^ 4 

pa = B^—gnn T^g • (25) 

The density of one neutrino species at that epoch is 

P.-^^T^NS, 

so that ^ 

6N, = 13.6g^lB^ = . (26) 

For the last equality we used ^^rh = 100, appropriate if Trh ~ 1 TeV. The alternative choice 
^(RH ~ 10 would reduce SNi, by a factor 2 or so. 

In models with more flatons, denoted by a label /, the quantity -Ba to be used in Eq. (|26|) 
is given by 

Ejnj (r(/ ^ aa) + ir(I ^ Xa)) 

E/^/r(/-.x) ' 

where rij is the number density of the flaton I just before reheating. 

If the flaton fields suffered negligible energy loss until they start their sinusoidal oscillation 
about their vev, one could in principle calculate the rij by solving the field equation of motion 
under the potential Eq. (|30|). The same thing is possible if the energy loss can be calculated, 
the only known paradigm for that purpose being the parametric resonance approximation. 
Such a calculation would be difficult, and its uncertainty impossible to quantify at present. 

One can however say something useful without knowing the n/, by considering the quan- 
tities 

If all of these quantities were equal, they would be equal to i?a- More usefully, if they are 
all known to have an upper or a lower bound r, in some regime of parameter space, then 
-Ba has the same bound. In the case of an upper bound, one concludes from Eq. (P^ ) that 
6N,y < 4.4r, making the model in this regime compatible with a given bound on 6N,j if r is 
small enough. In the case of an upper bound, one concludes that 6Ni, > 4.4r, making the 
model in this regime incompatible with a given bound if r is big enough. 
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3 Flaton and flatino spectrum 



3.1 The superpotential 

The model we consider contains two flaton superfields P and Q, interacting with the super- 
potential |Tl], [12 

W^flaton = -J^P^'^'Q ■ (29) 

We deal with the simplest case n = 1 and assign the PQ charges —1 and 3 to P and Q, 
respectively. Here we note that quantum gravity may break PQ symmetry, giving nonrenor- 
malizable terms which invalidate the PQ solution to the strong CP problem. A way to avoid 
this is to impose a certain discrete gauge symmetry forbidding sufficiently higher dimensional 
operators [0]. However, to have such a discrete symmetry, one has to extend the model 
beyond the simple superpotential (|29|) under consideration. Our analysis can be applied to 
such extended cases with a straightforward generalization. 

With the inclusion of the soft susy breaking terms and the cosmological constant, the 
potential is 



The soft parameters mp, uiq and Af are all of order 100 GeV in magnitude. It is assumed 
that mp and are both positive at the Planck scale. The interactions of 0p with the right 
handed neutrino superfields give radiative corrections which drive mp to a negative value 
at the PQ scale, generating vevs vp and vq for respectively |0p| and |0q|. According to 
Eqs. (|TD|) and (|n|), both vevs are roughly of order 10^° GeV. As we shall discuss in Section 
the radial oscillations of the flaton fields (f)p and (f)Q correspond to two flatons, while the 
angular oscillations correspond to a third flaton and the axion. 

3.2 Flaton spectrum 

We write the flaton fields as 

, Vp + P 

- ^e^-- (31) 

From now on, we shall take vp and vq as the independent parameters trading with m?p and 
tuq in the potential Eq. (^0|). The vevs are taken to real and positive, and we shall take the 
other independent paramerters to be and / to be real with opposite sign. 
The axion field is 

a = -^Ap + ?>^AQ (32) 

-rpQ rpQ 
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where FpQ = Vp + 9vq. The orthogonal field to the axion (both are CP Odd) corresponds 
to a fiaton. It is 

i.'=-^A^-S^A,. (33) 

-rpQ rpQ 

With our choice Aff<0, the vev is at ip' = 0, and the mass-squared is 

fAfVpF^^_ f 



where 



^ VpVQ 

9 ~ 2M, 



PI 



(35) 



X ^ ^. (36) 

For future convenience we have introduced a quantity /i, related to the g appearing only in 
the DFSZ model. At this stage results depend only on the ratio fi/g defined by (^) and 
they apply to both models. 

The other two fiatons correspond to the CP even fields P and Q. They have a 2®2 mass 
matrix whose components are 

M'qq = (37) 



„ „ vi Mix / f „ Ml 

Ml^x 9 + x2 V ^ 9 + ^ 

9 9wi(x2 + 3) Ml /•2 



Here two mass parameters mlp^rnq in Eq. (|30D are replaced in favor of vp,vq. Performing 
the rotation from the fiavor basis ||P QH to the mass basis F2II 

P = cos aF2 — sin aFi 

Q = sin + cos aFi , (38) 
we find the mixing angle a determined by 

Ml,p-Mi^Q_ x' + 3 



cos 2a 



M2 - m2 V9 + 42x2 + x^ 



'F2 ^''-'Fi 

sin 2a = — ^-f-^ = 39 
Ml^-M\ V9 + 42x2 + x4 ^ ^ 



where e = signMpg as we have the relation, M^^ - M\ = |M|,q/3x| V9 + 42x2 + a^. Later, 
the decay rates can be expressed in terms of cos 2a and sin 2a without ambiguity in fixing 
the angle a itself. The two eigenstates Fi, F2 have masses. 



Ml=^{ ^(12(x2 + 3)^ + (3 - x^)^) ± 1^(12^ + ^)|V9 + 42x2 + x4^ 
' ^ \9 9 1^ 9 9 1-^ J 



(40) 
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with Mf^ > Mf,. 

The requirement Mp_^ > gives the constraint 



where ^ 

yi,2 = (21 + x' ± V9 + 42x2 + 



5X- 

or 



^ (21 + x^ - V9 + 42x2 + x^) < -^-^ < ^ (21 + x^ + V9 + 42x2 + x^) . (42) 
One can find Mi^^ < M^/ for all the parameter space, and thus Fi is the lightest fiaton. 

3.3 Flatino spectrum 

Prom the superpotential VFflaton we can directly extract also the flatino 's mass matrix whose 
eigenvalues are 



Ml = ^^[x' + 2±2V^^] =94t/^2^x'^2±2A2^] (43) 
^2.1 4yx2^ 



The eigenstates Fi, F2 are related to the flavor states P, Q by 



F-i — cos aP + sin aQ 

F2 = - sin aP + cos (5(5 (44) 



where the angle a is determined by 

cos 2a — 



x^ 



Vl + a: 

A parameter space analysis indicates that we have always Mf-^ < 2Mp_^. This automatically 
forbids the decay of Pi to flatinos leaving open only the decay into flatinos of the heavier P2 
and ■0' flatons. 

4 Decays involving only flatons, flatinos and axions 

In this section, we analyze the various decay rates between flatonic fields. We begin with 
the decay channels induced by the kinetic term and the superpotential H^flaton, which are 
common to the KSVZ and DFSZ models. 
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4.1 Derivative and cubic interaction terms between flatons 

The flaton interaction terms with at least one derivative are given by the Lagrangian 



' ''"2^ (9af + If mf - e'-^a^'da] (46) 

^^Q -^PQ -^PQ -^PQ 

From this expression we can extract the following terms expressed in mass eigenstates: 

(i) the trilinear derivative interactions with no axions, 

(dib')'^ [f-9 sin a + cos a) Fi + (q cos a + x^ sin a) F2I ; (47) 

FpQxVx"^ + 9 ^ ^ ^ ^ 

(ii) the trilinear derivative interactions with only one axion, 

Lp-ip'a = dijj'da [(cos a — x sin a) F2 — (sin a + x cos a) Fi] ; (48) 

-^pQVa;^ + 9 

(iii) the trilinear derivative interactions with two axions, 

Lp aa — \dixa\'^ ((9 cos g — X sin a) Fi + (9 sin a + x cos a) F2) . (49) 

-^PQ V ^2 + 9 

All the above derivative interactions can be transformed in scalar interactions if we are 
working at tree level and with on-shell external particles 

4>i {dM (d^h) = \ {Ml - Ml - Ml) 010203 (50) 
The cubic interactions come also from the superpotential and the soft terms 

Ml, ^ 2^ Ml, vlix^ + 9)) '^[2^ Ml, vp {x^ + 9) J ^ 



PQ. ( 3/ p^y + i LM:^\ g^y . (51) 



In the mass basis, we get 



r Ml,, siiio- / 2 2 \ 

JF2-Fii^i = 1 3 ' . : (—2a; COS a + a; sin a + cos a sin a 1 + 

[ FpQxV9 + x2 ^ ^ 

6 p u^^/x^ + 9 / -,„2- ^ ■ ■\ 0S2 

( — 18 cos a sm q; a; + 9 sm o; a; + 3 cos a x 

g^xFpQ ^ 

- cos a sm' a{-9 + x')) } F^F^ = ^^^F^F^ (52) 
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For the full trilinear F2ip'ip' interaction, we have to add up the terms in Eqs. ( |17| ) and ( pT| ) 
to obtain 



2xV9 + x^FpQ 
(9 cos a + sin ax^) f 1 — 2— | | 



^— (3 sin a + X cos a) (^x^ + 9^ 



(53) 



Collecting the above formulae one finds the decay rates among flatons and axions 



r (F2 — > aa) = 
r (Fi — > aa) = 

r (F2 ^ FiFi) 
r (F2 ^ v-V) 

r (F2 ^ a7/>') = 
r (V-' ^ aF2) = 
r(V>'-.aFi) = 



327rF2Q (x2 + 9) 

327rFp2Q (x2 + 9) 
1 



77^ (xcosa + 9sma) 



-X sin a + 9cosa)^ 



1 



12 

1 — 4^T2- |^F2FiFi I 



1-4- 



1 Ml ( Ml,\\ 

'1 — I (3 COS a — 3x sin a) 



167rFp2Q(x2 + 9) V 



1 



2 \3 



167rF2Q(x2 + 9) V Ml,) 

r2 \ 3 



(3 cos a — 3xsina) 



Ml. 



„ , ^ -(1 ^] (3 sin a + 3a; cos 



(54) 
(55) 

(56) 

(57) 
(58) 
(59) 
(60) 



Energy conservation will of course forbid some of these reactions, depending on the flaton 
masses. As Mp^ < M^' the channels Fi —>■ tp'ip' and Fi —>■ ip'a are always forbidden. 



4.2 Interaction terms between flatons and flatinos 

The trilinear Lagrangian terms responsible for the decay of flatons or flatinos are 

-iMpi \ _^PQ i^PQ J 

'2PvpPQ + 2%^ (3t;Q^' + vpa) P^^Q] 



2M, 



(61) 



PI 



PQ 



(the tilded fields are the fermionic superpartner of the respective P and Q scalars). Let us 
denote the Yukawa couplings between the flaton (or the axion) and the flatinos in mass basis 
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by —Lyuk = Yijk'PiFj (1; lb) Fkl'^ wheie 75 is taken for 0j = a, ifj' . We find from Eq. (|6T|) the 
following expressions for the Yukawa couplings 

6//iVa;2 + 9 • x 2- • • o~i 

r P P s = (xcosa — sma) cos a — x sm a sm 2a 

6iW^^T9 . s . 2 ~ , • • o~i 

p.Po = (xcosa — sma) sm a + xsmasm/a 

GffiVx^ + 9 . 1.. ^ • o~i 

r P p P = — a;sinacos2a H — (sma; — xcosa) sm/a 

Gffiy/x^ + 9 . N 2 ~ , • o~l 

r p p p = (xsma + cosa) cos a + a; cos a sm 2a 

6//iVx2 + 9 . ^ . 2 ~ -0-1 

^F2F2F2 ~ 'p sm a + cos a j sm a — x cos a sm 2aJ 

Qfa^/x"^ + 9 r „^ 1 . 
r P P P = +a;cosacos 2a (cosa + xsma) sm 2a 

Y ,,p p = [— (3 + x^) cos^ a — 3xsin2a] 

^V'i^2F2 = ~^^F~^~ ^ ^ ^"^^ ^^"^^ a + 3x sin 2a] 
^V'i^i^2 = ^^F^ [-3xcos2a + ^ (3 + x^) sin 2a] 
K p p = ^ f2a; cos^ a — x^ sin 2a] 
p = [2x sin^ a + sin 2a] 

^aFiFa = ^{.^ [-^"^ COS 2a - X siu 2a] (62) 
fl'a^-TpQ 

From this we can extract the decay rates for Fi — > FjFi^, ip' —>■ FjFf^, or F2 FiFi [tp', a) 

^(--^^-) ^ 



16n \ Ml \ Ml 

2 /i#2 
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1 1 

X 

r(F2^Fia) = 

where 5* is a symmetric factor (1/2 for identical final states or otherwise 1). 

5 Interaction of flatons and flatinos with matter fields 

Now we study the interactions of the flatons with matter and supermatter. Through these 
interactions the flatons and flatinos decay into ordinary matter and axions, and the produc- 
tion of the latter must be sufficiently suppressed to satisfy the nucleosynthesis limit on the 
effective number SN^, of extra neutrino species. (In this paper, we are for simplicity assuming 
that no flatino is stable.) 

5.1 KSVZ model: Interactions between flatons and gluons 

In the KSVZ (hadronic) model, the interaction with matter is 

W^flaton-matter = flEiEiE'^P (68) 

where Ei and Ef are additional heavy quark and antiquark superfields. 

The only decay mode available for the flatons is into two gluons coming from the anomaly 
(when the space phase will be available, we have to take into account also the decay into 
massive gluinos, in this discussion we neglect such a possibility). The respective one loop 
corrected decay rates are 

where Ne is the total number of the superheavy exotic quark fields {Me = fiEVp ^ Mp-). 
We do not consider the flatino decay into a gluon and a gluino which will be irrelevant for 
our discussion. 




(66) 
(67) 
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5.2 DFSZ model: Interactions between flatons/flatinos and ordi- 
nary matter 

In the DFSZ model, the interaction is 



W^flaton-matter = ^NNP + -^H,H,pQ 



(70) 



where are the right handed neutrino superfields and i/1.2 the two Higgs doublets. Due to 
the second term we can provide a solution to the \x problem |^ . In such case we can add to 



the superpotential of the minimal supersymmetric standard model also the terms hulH2N 
that generate the necessary mixing between left and right neutrinos to implement a see-saw 
mechanism which can explain the solar and atmospheric neutrino deficits. 

The decay properties of the flatons now involve the direct interactions between flatons and 
ordinary matter and supermatter In general the interaction between flatons and Higgs fields 
are quite interesting due to the fact that these two sectors, after the spontaneous breaking 
of the PQ and the EW symmetry, mix together. We notice that the Peccei-Quinn symmetry 
prevents the introduction of a SUSY invariant mass term fiHiH2, solving automatically the 
so-called n mass problem as mentioned before. 

Let us start by writing the Higgs-flaton potential 



ViH, 



+ gHiH2 Ag 



Ml 



PI 



M, 



+ 3/ 



PI 




+ c.c. 



(71) 



When the fields (/>p,q get vevs, the m\HiH2 mass term is generated dynamically. The size 
of such a term is fixed by 



mg = /i {Ag + [x^ + 3^ 



(72) 



In the limit \mW ^ the masses of the pseudoscalar of the CP even scalar Higgs field 
and of the charged Higgs fields are almost degenerate 



rrin 



sin P cos P 

so from the constraint of positivity of such a masses we get 

9 



(73) 



(74) 
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In such a limit we also know that the mass eigenstate of the CP even electroweak sector 
H^, hP and of the CP odd one A^, are 



= - sin (3 hi + cos (3 hi 

h^ = cos/3/i? + sin/5/i° 

A° = sin/5 A° + cos/3 A° 

G° = cos (3 - sin (3 



(75) 



where Hi = {vi + h^ + iA^) and H2 = ^ (f 2 + + ^^2) the gauge eigenstates and 
tan (3 = V2/V1. To allow the flaton decay into A^, we want it to be light so that small tan/3 
is preferred in our discussion. Hereafter we will take tan (3 = 1. 
From Eq. (|72|), we find 



V^,, = ^/i^ {hf + hf + Af + Af' ^ ^ 



Ag^I 



+C.C. 



2' 
'P_ 

Vp 



Q_ 



PQ 



9 



f 2(P , Q 



Vp 



p 



i)' + x^-^i^ 4 h i 



PQ 



Vp 



PQ 



(76) 



It is then simple manner to get the decay rates for the kinematically more favorable decay 
channels -Fi 2 ^ h^h^ and ip — > h^A^ 



r(F, 
r(F, 
r 



(x^ + 9) /x^ 
327rF|Q 16x2 M|,^ 



1 - 



4M|o' 



1/2 



F^_hh\ 



327rF|Q 16x2 



1 - 



AM, 



1/2 



fx- 



F2 



\A 



F2hh\ 



(77) 



1 - 



(Mfeo - M^o) 



2\ 1/2 



1 - 



(Mfeo + M^o) 



2\ 1/2 



I A 



where 



Apihh 
Apihh 



sin2^[ 

sin2/3[ 
'A 



X cos a — sin a) — 4x^ — sin a] + 2 (x cos a — sin a) 

.2f 



(x sin a + cos a) + 4x — cos a] + 2 (x sin a + cos a) 

9 



A^'hA = I — - 6- 



X 



(7J 



The flatino decay into ordinary particles comes from the superpotential W = ■j^HiH2pQ. 
We find that the flatino decay into a Higgs and a Higgsino (more precisely, the lightest neu- 
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tralino Xi) has the rate; 

(79) 

where Cp^ — (sina + x ^cos(5)A^^^ and Cp^ — {cosa — x ^ sm.a)N-^^. Here N-^^^ denotes the 
fraction of hghtest neutrahno in Higgsinos. 

Let us now consider the flaton decay into ordinary fermions or sfermions. The mixing 
terms between flaton and Higgs fields allow a direct tree level coupling (after full mass matrix 
diagonalization) between the usual fermions and flatons. Parameterizing such a mixing with 
the parameter 9pH the effective flaton-fermion interaction is hfOpH so that the rate of decay 
is 

r(F.^/ + /)=7V.^M..(l-4^j (80) 
where Nc is a color factor for the fermion /. Since 6fh — (^^) 

r (f, ^ / + /) /r (F, ^ a + a) ~ h)vlw/Ml ~ mJ/M|^ < 1 . (81) 

Therefore, the rate of the flaton decay into ordinary fermions cannot be made sufficiently 
larger than that into axions. 

For the coupling between sfermions and flatons, we have two contributions. One is a 
direct coupling coming from the scalar potential 



= h [hdt^^PDL Dr + hetan^EL Er + KUl Ur ) 82 

FpQ a; V2 ^ ^ 

(Fi {x cos a — sin a) + F2 (cos a + a; sin a) ) + h.c. 

where D* denote down-type squark, etc. 

The other arises from an indirect coupling induced by the mixing between Higgs and 
flaton fields as for the fermion case. Taking in consideration the cubic soft A-terms we find 

K/ / = hdAd Of, h, F^DlDr + h.A, Op, h, F^ElEr + KA^ Op, h, FJJlUr (83) 

so that effectively we have couplings of the size 

Gpsfermion hf{lJ, + Af)^^ (84) 



Diagonalizing the sfcrmion mass matrix we can write /r/l = oii/i/i + ^22/2/2 + 012/1/2 
(where an oc hf so for hf we have 012 1). Considering the decay of the light flaton 
we get 

r(F.^/. + /,).iV.^aiJl-4^ (85) 
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Table 1: Direct decay channels involving only fiatons (/ = Fi, F2 and ip'), flatinos [Fi and 
F2), and the axion (a). The only decays which must occur are Fi aa, ip' — >■ aFi, and 
either ip' —>■ aF2 or F2 — > tp'a. Any of the others may be forbidden by energy conservation, 
if the decaying particle is too light. 



Fi F2 




Fi F2 


aa aa 


aFi 


none Fil 




aPs 






p.p. 




FiF^ 






ip'ip' 







As observed in Ref. |jT5|, the fiaton may decay efficiently to two light stops as /ij ~ 1 and 
aij ~ 1 and thus a large splitting between light and heavy stops helps increasing the fiaton 
decay rate to light stops. This kind of mass splitting occurs also in the Higgs sector and 
furthermore the light Higgs {h°) is usually substantially lighter than the heavy Higgs {H°) in 
the minimal supersymmetric standard model. This should be contrasted to the case with the 
mass splitting for stops which requires some adjustment in soft parameters. In this paper we 
concentrate on the fiaton decay into Higgses, as the probably dominant mode, which in any 
case provides a lower bound on the decay rates to ordinary matter and therefore an upper 
bound on 5N,,. 



6 Parameter space analysis 

We have now evaluated all of the direct decay rates for fiatons and flatinos into channels 
involving axions, as summarized in Table 1. We have also evaluated some of the contributions 
to the direct decay rates of fiatons and flatinos into hadronic matter X. The ultimate 
objective is to evaluate the decay rates r(/ X), r(/ Xa), and r(/ aa) for each of 
the three fiaton species, so as to evaluate SN^ through Eqs. (|^) and (|^). For the reactions 
Fi — >■ aa and Fi — >■ X, the direct rates are the same as the total rates, but in the other cases 
one has to consider also chain reactions. The reactions F2 — > X and X can go either 

directly or through chains, and their rates are 

r(^'^x) = rdir(V''^x) + r(v^'^FiFi) 

+r(V'' FiF2)B{F2 ^ X) + r(V^' ^ F2F2)B\F2 ^ X) 

r(F2^x) = rdir(F2^x) + r(F2-^FiFi) 

+r(F2 -> F,F2)B{F2 ^ X) + r(F2 ^ F2F2)B\F2 ^ X) , (86) 
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where B denotes a branching ratio. The reactions ip' ~^ aX, F2 aX can go only through 
chains, and their branching ratios are 

T{tP' ^ aX) = T{ij' aF2)B(F2 ^ X) + r(^' FiF2)B(F2 ^ Fia) 

+2T{ij' ^ F2F2)B{F2 ^ Fia)B(F2 ^ X) (87) 

r(F2 ^ aX) = r(F2 a^')B(V^' ^ X) + r(F2 FiF2)B(F2 ^ Fia) 
+21(^2 F2F2)B{F2 Fia)B(F2 X) 
+2r(F2 ^ ^'i)')B{i)' aX)B{^lj' X) . 

The rates for producing more axions are hkely to be suppressed, because every term in the 
analogous expressions containing the product of more than two or more branching ratios. 



Through Eqs. (^q) and (^TQ , these expressions provide the basis for a calculation of 6N,y, 
as a function of the parameters, if the relative values of the initial flaton densities rij are 
known. Here, we perform the less ambitious task, of trying to identify a region of parameter 
space allowed by an upper bound SN^, of order 1 to 0.1. (Recall that the former bound is 
roughly the present one [Q, while improvements in the foreseeable future will either tighten 
the limit by an order of magnitude, or else detect a nonzero 6N^.) We shall find that such 
a region probably does not exist in the KSVZ case, but that it certainly does exist in the 
DFSZ case. 

6.1 Parameter space of KSVZ models 

As discussed already, the flatons Fi 2 can decay into two gluons, so let us try to see if the 
generous requirements -8^1,2 = r(-Fi 2 — > aa)/r(Fi 2 — * gg) < 3 [see Eq. (^)] can be fulfilled 
in any region of parameter space. From Eqs. (^41), ( ^Sf ) and (|69|), we get 

Bp^ ~ 8 X 10-^iV| 
B^^ ~ 8 X 10-^Ar| 



sin a (x^ + 9) 
X (9 cos a — X sin a) 

cos a (x^ + 9) 
X (9 sin a + xcosa) 



for as (Mp) ^ 0.1. With e = +1, we find that Bp^ < 7 x 10 "^N^ for any value of x, and 
Bp^ < 6 X lO-^A^I/x^ in the limit x ^ 0. Therefore, Bp^ can be made small enough for 
X ~ 0.1 but Bp-^ cannot taking a reasonable value of Np- A way out would come from a 
cosmological evolution of the flatons. That is, one could imagine a situation in which the 
flatons oscillate only along the direction of F2 so that the populations of Fi and ip' (which 
can decay into aFi) are suppressed by the order of 10^ compared to that of F2. But this is 
not probable. On can find the similar behavior for e = —1 in which case Bp^ can be made 
small in the limit x — > 0. 

6.2 Parameter space of DFSZ models 

In this subsection, we will try to find a region of parameter space of the DFSZ models. For 
this, we consider rather stringent requirement; Bj < 0.1 for all three flatons and two fiatinos. 
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According to Section p75t this requirement will ensure that dNi, < 0.1. Our strategy is to 
first write down a set of conditions which ensure the requirement, and then to identify a 
region of parameter space in which these conditions are satisfied. 



6.2.1 The conditions 



We first note that the decay rates calculated in the previous sections are functions of the 
4 variables x = vp/vq, f /g, Aj/fi and Ag/fi disregarding their overall dependence on Fpq. 
To be as independent as possible of the soft supersymmetry breaking parameters we will 
try to make analytic computations on the rates of the flaton decays into Higgs particles, in 
particular into the lightest Higgs boson (h^) whose mass has an automatic upper bound of 

< 1 and X > 1. 



140 GeV [Q. We will concentrate on the region with ^ negative and 
To open the decay channels of the flatons into Higgs particles we have, in particular, to 



impose M^/ > Ma > which requires — < 



X 



with 



9 



x^ >2 



The positivity of flaton masses requires 



9 , f 2 
- H — X 



(89) 







f_ 


< x^ 








9 



(90) 



In case the fiatino production rates are sizable, we also impose Rp^ = T{F2 Xih^)/^{^2 
aFi) > 10 and Mp^ > M^j + M^o to open the decay mode Fi Xih^- These two conditions 
give rise to the restrictions; 



|//^| < 0.02xA^xi and 



g 6x 



(91) 



For the latter condition, we required Mp^ > /i. 

Then we study, in our limit, the constraints given by the conditions R^i = F (■?/'' 
/iM°)/F(V'' ^ aFi) > 10 and Rp, = T{Fi h^h^)/V{Fi aa) > 10. The ratios R are 



Rih' 



R 



R 



■F2 



1 



' \ 2 2 

9\ 



144 V/j ^1 



4M4//i g 
,4 



9 o / 3,2 ^ 2)2 



(92) 



10" 



(^ + 18- +2) 
9 



where Mp^ ~ 

121^1. 

' 9 ' 



9 M 



2 

F2 



12 4a;V for ^ < 121^1, and M, 



F, ^ for ^ > 
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6.2.2 A viable region of parameter space 



Now we identify a region of parameter space, in whicli tlie above conditions are all satisfied. 
Recall that we are considering only the region ^ 1, f /g negative, and \ f / g\ <C 1. Within 
this region, we consider the four regions 



I) 
II) 

III) 
IV) 



A, 


> 


f_ 










9 



X 



2 < 

/ 
9 
A 



A, 


< 


f_ 










9 



X 



< 



< 



< 2 



Depending on ^ < 12|^| or ^ > 12|^| we define regions a or region b. 

We find that all of the a regions are forbidden, and so is the IVi, region. The constraints 
for the other regions are as follows. 



X > 14, Ag < 0, 12 
A 



I 


< 






9 







< 



1^ 

x>9, Ag>0, 

A. 

2 < 



x'^ 
— < 
2 



A„ 



< 2 



(93) 



< 3 X 10 



-2 







/ 




2 < 


^9 




/ 


< 




^^ 




< 








9 










9 



X 



(94) 



Ilh Ag>0, 1< 2 



X 



< 



At 



< 2 X 10" 



/ 


< 








9 







< 2. 



(95) 



In addition to these conditions, we need to add the original requirements, that f/g he 
negative with \ f /g\ <^ 1. (The latter condition is automatically satisfied in regions lib and 
I lib, but not in region /;,.) The other original requirement a; ^ 1 is automatically satisfied 
in all three regions. Note that all cases one requires x'^\f/g\ ^1. 

Within each of these allowed regions, the parameters can take on their natural values 
1^4/1 ~ \Ag\ ~ ~ 100 GeV, vp ~ vq and |/| ~ \g\ ~ 1, within a factor 10 or so. 
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7 Conclusions 



We have explored the cosmology of a supersymmetric extension of the Standard Model, 
which has a Peccei-Quinn symmetry broken only by two 'fiaton' fields 0p and 0q. They 
have two radial modes of oscillation and one angular mode (plus the axion), corresponding 
to three fiatons with mass of order 100 GeV, and there are two fiatinos with roughly the 
same mass. The fiatons are the generalizations of the saxion which appears in non-flaton 
models, and the fiatinos are the generalizations of the axino. 

The fiaton models have the virtue that FpQ is predicted in terms of the electroweak 
scale and the Planck or GUT scale. In the canonical model we have estimated FpQ ~ 
-|^gio.4±o.9 QgY^ g^]2d in the more complicated ones Fpq > 10^^'^ GeV. 

We have assumed that (pp has a positive effective mass-squared in the early Universe, so 
that thermal inflation occurs, allowing rather definite predictions. The axion is a good dark 
matter candidate in all cases. (By good dark matter candidate we mean one whose density is 
predicted to be roughly in the right ballpark.) Bearing in mind that PQ strings arc produced 
after thermal inflation, this simply corresponds to the received wisdom in the canonical case 
where FpQ ~ 10^° GeV. But the axion is also a good candidate when Fpq is bigger, because 
entropy production from fiaton decay more than compensates for the increased axion density 
before fiaton decay. Besides the axion, the LSP is also a good dark matter candidate in the 
canonical model. (In the other models, the reheat temperature is too low to thermalize the 
LSP, which must therefore be unstable.) In KSVZ models, a third good candidate is the 
supermassive particle whose mass comes mostly from a coupling to one of the fiatons. 

Our main concern has been with the highly relativistic axion population that is produced 
by fiaton decay It remains relativistic to the present and therefore makes no contribution 
to the dark matter, but it is dangerous for nucleosynthesis because it is equivalent to very 
roughly 6Ni, ~ 1 extra neutrino species. At present, the bound at something like the 2-cr 
level is 6Ni, < 1.8 for the 'high' deuterium nucleosynthesis scenario, and 6Ni, < 0.3 for the 
perhaps favored 'low' deuterium scenario. In the foreseeable future one will either have a 
bound SN,^ ^0.1, or a detection of SN^,. 

We have calculated the rates for all relevant channels and examined the constraint that 
the energy density of these axions does not upset the predictions of the standard nucleosyn- 
thesis. We confirm the earher conjecture, that the KSVZ case is probably ruled out even 
by the present bound ^A^^ ^ 1. For the DFSZ case there are more decay channels. To 
evade complicated phase space suppressions we concentrate on the decay of the fiatons into 
Higgses, as the mass of the lightest Higgs boson has naturally a relatively low upper bound 
and the mass of the other Higgs boson and fiaton fields are fixed by the parameters of the 
flatonic potential itself. In this way, we have found a region of parameter space, in which 
5Ni, will certainly be ^0.1, and in which the parameters can take on their natural values 
within a factor ten or so. In its DFSZ version, the fiaton model of PQ symmetry breaking 
will be a candidate for explaining a future detection of nonzero SN^,. 

An interesting question, lying beyond the present investigation, is whether the allowed 
region of parameter space can be achieved in a super gravity model with universal soft pa- 
rameters. 
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